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What is claimed is: 

1 . A computer readable medium or media, comprising: 

(a) a data structure relating a plurality of reactants to a plurality of reactions of 
a biochemical reaction network, 

wherein each of said reactions comprises a reactant identified as a substrate of 
the reaction, a reactant identified as a product of the reaction and a stoichiometric 
coefficient relating said substrate and said product, and wherein at least one of said 
reactions is a regulated reaction; and 

(b) a constraint set for said plurality of reactions, 

wherein said constraint set comprises a variable constraint for said 
regulated reaction. 

2. The computer readable medium or media of claim 1, wherein said 
variable constraint is dependent upon the outcome of at least one reaction in said data 
structure. 

3. The computer readable medium or media of claim 1, wherein said 
variable constraint is dependent upon the outcome of a regulatory event. 

4. The computer readable medium or media of claim 1, wherein said 
variable constraint is dependent upon time. 

5. The computer readable medium or media of claim 1, wherein said 
variable constraint is dependent upon the presence of a biochemical reaction network 
participant. 

6. The computer readable medium or media of claim 5, wherein said 
participant is selected from the group consisting of a substrate, product, reaction, 
protein, macromolecule, enzyme and gene. 
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7. The computer readable medium or media of claim 1, wherein said 
biochemical reaction network comprises metabolic reactions. 

8. The computer readable medium or media of claim I, further 
comprising a regulatory data structure, wherein said variable constraint is dependent 
upon an outcome of a regulatory event represented by said regulatory data structure. 

9. The computer readable medium or media of claim 8, wherein said 
regulatory data structure represents a regulatory event selected from the group 
consisting of transcription of a gene, translation of an RNA, post-translational 
modification of a protein, inhibition of a protein, activation of a protein, assembly of a 
protein, change in pH, change in redox potential, change in temperature, passage of 
time, and degradation of a protein. 

10. The computer readable medium or media of claim 8, wherein said 
regulatory event is due to a signal transduction pathway. 

1 1 . The computer readable medium or media of claim 8, wherein said 
biochemical reaction network and said regulatory data structure represent reactions or 
events that occur in a single cell. 

12. The computer readable medium or media of claim 8, wherein said 
biochemical reaction network represents reactions that occur in a first cell in a 
population of cells and said regulatory data structure represents events that occur in a 
second cell in said population. 

13. The computer readable medium or media of claim 12, wherein 
said population of cells comprises cells of a multicellular organism. 

14. The computer readable medium or media of claim 1, further 
comprising a constraint function that correlates an outcome of a regulatory event with 
said variable constraint. 
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15. The computer readable medium or media of claim 14, wherein 
said constraint function is binary. 

1 6. The computer readable medium or media of claim 14, wherein 
said regulatory event is represented by Boolean logic. 

1 7. The computer readable medium or media of claim 1 , further 

comprising 

(c) commands for determining at least one flux distribution that minimizes or 
maximizes an objective function when said constraint set is applied to said data 
structure, wherein said at least one flux distribution determines a systemic property of 
said biochemical reaction network, wherein said systemic property is dependent upon 
the flux through said regulated reaction. 

18. The computer readable medium or media of claim 17, wherein 
said commands determine a range of feasible flux distributions that minimize or 
maximize an objective function when said constraint set is applied to said data 
representation. 

19. The computer readable medium or media of claim 17, wherein 
said commands comprise an optimization problem. 

20. The computer readable medium or media of claim 19, wherein 
said optimization problem comprises a linear optimization problem or a nonlinear 
optimization problem. 

2 1 . The computer readable medium or media of claim 1 7, further 
comprising a user interface capable of sending at least one command for modifying 
said data structure, said constraint set or said commands for applying said constraint 
set to said data representation, or a combination thereof. 
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22. The computer readable medium or media of claim 2 1 , wherein 
said user interface further comprises links which a user may select to access 
additional information relating to said plurality of reactions. 

23. The computer readable medium or media of claim 1, wherein said 
data structure comprises a set of linear algebraic equations. 

24. The computer readable medium or media of claim 1 , wherein said 
data structure comprises a matrix. 

25. The computer readable medium or media of claim 1, further 
comprising commands for representing said at least one flux distribution as a flux 
distribution map. 

26. The computer readable medium or media of claim 1, wherein at 
least one reactant in said plurality of reactants or at least one reaction in said plurality 
of reactions is annotated. 



27. The computer readable medium or media of claim 26, wherein 
said annotation comprises assignment of said at least one reactant to a compartment. 

28. The computer readable medium or media of claim 27, wherein a 
first substrate or product in said plurality of reactions is assigned to a first 
compartment and a second substrate or product in said plurality of reactions is 
assigned to a second compartment. 

29. The computer readable medium or media of claim 26, wherein 
said annotation comprises assignment to an open reading frame or protein. 

30. The computer readable medium or media of claim 26, wherein 
said annotation comprises a confidence rating. 
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3 1 . The computer readable medium or media of claim 1 , further 
comprising a gene database relating one or more reactions in said data structure with 
one or more genes or proteins in particular organism. 

32. The computer readable medium or media of claim 1, wherein said 
biochemical reaction network comprises reactions that are selected from the group 
consisting of glycolysis, the TCA cycle, the pentose phosphate pathway, respiration, 
biosynthesis of an amino acid, degradation of an amino acid, biosynthesis of a purine, 
biosynthesis of a pyrimidine, biosynthesis of a lipid, metabolism of a fatty acid, 
biosynthesis of a cofactor, metabolism of a cell wall component, transport of a 
metabolite and metabolism of carbon, nitrogen, sulfur, phosphate, hydrogen or 
oxygen. 

33. The computer readable medium or media of claim 1, wherein a 
plurality of said reactions are regulated reactions and said constraints for said 
regulated reactions comprise variable constraints. 

34. A method for determining a systemic property of a biochemical 
reaction network, comprising: 

(a) providing a data structure relating a plurality of reactants to a plurality of 
reactions of a biochemical reaction network, 

wherein each of said reactions comprises a reactant identified as a substrate of 
the reaction, a reactant identified as a product of the reaction and a stoichiometric 
coefficient relating said substrate and said product, and wherein at least one of said 
reactions is a regulated reaction; 

(b) providing a constraint set for said plurality of reactions, 

wherein said constraint set comprises a variable constraint for said regulated 
reaction; 

(c) providing a condition-dependent value to said variable constraint; 

(d) providing an objective function, and 

(e) determining at least one flux distribution that minimizes or maximizes said 
objective function when said constraint set is applied to said data structure, 
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thereby determining a systemic property of said biochemical reaction network. 

35. The method of claim 34, wherein said value provided to said 
variable constraint changes in response to the outcome of at least one reaction in said 
data structure. 

36. The method of claim 34, wherein said value provided to said 
variable constraint changes in response to the outcome of a regulatory event. 

37. The method of claim 34, wherein said value provided to said 
variable constraint changes in response to time. 

38. The method of claim 34, wherein said value provided to said 
variable constraint changes in response to the presence of a biochemical reaction 
network participant. 

39. The method of claim 38, wherein said participant is selected from 
the group consisting of a substrate, product, reaction, enzyme, protein, macromolecule 
and gene. 

40. The method of claim 34, wherein said biochemical reaction 
network comprises metabolic reactions. 

41. The method of claim 34, further comprising a regulatory data 
structure, wherein said value provided to said variable constraint is changed due to an 
outcome of a regulatory event represented by said regulatory data structure. 

42. The method of claim 41, wherein said regulatory event is selected 
from the group consisting of transcription of a gene, translation of an RNA, post- 
translational modification of a protein, inhibition of a protein, activation of a protein, 
assembly of a protein, change in pH, change in redox potential, change in 
temperature, passage of time, and degradation of a protein. 
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43. The method of claim 41, wherein said regulatory event is due to a 
signal transduction pathway. 

44. The method of claim 41, wherein said biochemical reaction 
network and said regulatory data structure represent reactions or events that occur in a 
single cell. 

45. The method of claim 41, wherein said regulatory event comprises 
a regulatory reaction. 

46. The method of claim 41, wherein said biochemical reaction 
network represents reactions that occur in a first cell in a population of cells and said 
regulatory data structure represents events that occur in a second cell in said 
population. 

47. The method of claim 46, wherein said population of cells 
comprises cells of a multicellular organism. 

48. The method of claim 41, further comprising a constraint function 
that correlates an outcome of a regulatory event with said variable constraint. 

49. The method of claim 48, wherein said constraint function is 

binary. 

50. The method of claim 48, wherein said regulatory event is 
represented by Boolean logic. 

5 1 . The method of claim 48, wherein said constraint function 
correlates a first set of outcomes of said regulatory data structure with a first binary 
value and a second set of outcomes of said regulatory data structure with a second 
binary value. 
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52. The method of claim 48, wherein said constraint function 
correlates a set of outcomes of said regulatory data structure with a single integer 
value. 

53. The method of claim 34, wherein said flux distribution is 
determined by optimization. 

54. The method of claim 53, wherein said optimization comprises 
linear optimization or non linear optimization. 

55. The method of claim 34, further comprising a step of modifying 
said data structure or said constraint set, or both. 

56. The method of claim 34, wherein said data structure comprises a 
set of linear algebraic equations. 

57. The method of claim 34, wherein said data structure comprises a 

matrix. 

58. The method of claim 34, further comprising a step of producing a 
flux distribution map. 

59. The method of claim 34, wherein said biochemical reaction 
network comprises reactions that are selected from the group consisting of glycolysis, 
the TCA cycle, pentose phosphate pathway, respiration, biosynthesis of an amino 
acid, degradation of an amino acid, biosynthesis of a purine, biosynthesis of a 
pyrimidine, biosynthesis of a lipid, metabolism of a fatty acid, biosynthesis of a 
cofactor, metabolism of a cell wall component, transport of a metabolite and 
metabolism of a carbon source, nitrogen source, oxygen source, phosphate source, 
hydrogen source or sulfur source. 
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60. The method of claim 34, wherein said systemic property is 
selected from the group consisting of growth, energy production, redox equivalent 
production, biomass production, production of biomass precursors, production of a 
protein, production of an amino acid, production of a purine, production of a 
pyrimidine, production of a lipid, production of a fatty acid, production of a cofactor, 
production of a cell wall component, transport of a metabolite, development, 
intercellular signaling, and consumption of carbon nitrogen, sulfur, phosphate, 
hydrogen or oxygen. 

61. The method of claim 34, wherein said systemic property is 
selected from the group consisting of degradation of a protein, degradation of an 
amino acid, degradation of a purine, degradation of a pyrimidine, degradation of a 
lipid, degradation of a fatty acid, degradation of a cofactor and degradation of a cell 
wall component. 

62. The method of claim 34, wherein said variable constraint 
comprises a condition-dependent constraint value and a constraint function, wherein 
said variable constraint is modified by said constraint function acting upon said 
condition-dependent constraint value. 

63. The method of claim 62, wherein said constraint function is 

binary. 

64. The method of claim 34, further comprising providing a gene 
database relating one or more reactions in said data structure with one or more open 
reading frames or proteins in a particular organism. 

65. The method of claim 64, further comprising identifying an open 
reading frame that encodes a protein that performs a reaction in said plurality of 
reactions. 
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66. The method of claim 64, further comprising identifying a protein 
that performs a reaction in said plurality of reactions. 

67. A method for determining a phenotype of a mutant of an 
organism, comprising: 

(i) identifying a reaction that is not naturally present in a particular organism, 

and 

(ii) determining a systemic property of a biochemical reaction network 
according to the method of claim 34, wherein said data structure relates a plurality of 
reactants for said organism to a plurality of reactions of a biochemical reaction 
network of said organism and further comprises said reaction that is not naturally 
present in said organism. 

68. A method for determining a phenotype of a mutant of an 
organism, comprising: 

(i) identifying a reaction that is related to an open reading frame or protein in 
said gene database, and 

(ii) determining a systemic property of a biochemical reaction network 
according to the method of claim 34, wherein said reaction that is related to said open 
reading frame or protein is not present in said data structure or is constrained to have 
no flux. 

69. A method for determining the effect of a drug on the activity of 
one or more reactions in a biochemical reaction network, comprising: 

(i) identifying a reaction that is related to an open reading frame or protein in 
said gene database; 

(ii) identifying a candidate drug that alters expression of said open reading 
frame or activity of said protein, and 

(iii) determining a systemic property of a biochemical reaction network 
according to the method of claim 34, wherein said reaction that is related to said open 
reading frame or protein is not present in said data structure, is constrained to have a 
reduced flux, or is constrained to have no flux. 
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70. The method of claim 34, wherein a plurality of said reactions are 
regulated reactions and said constraints for said regulated reactions comprise variable 
boundary values. 

71 . A method for determining a systemic property of a biochemical 
reaction network at a first and second time, comprising: 

(a) providing a data structure relating a plurality of reactants to a plurality of 
reactions of a biochemical reaction network, 

wherein each of said reactions comprises a reactant identified as a substrate of 
the reaction, a reactant identified as a product of the reaction and a stoichiometric 
coefficient relating said substrate and said product, and wherein at least one of said 
reactions is a regulated reaction; 

(b) providing a constraint set for said plurality of reactions, wherein said 
constraint set comprises a variable constraint for said regulated reaction; 

(c) providing a condition-dependent value to said variable constraint; 

(d) providing an objective function; 

(e) determining at least one flux distribution at a first time that minimizes or 
maximizes said objective function when said constraint set is applied to said data 
structure, 

thereby determining a systemic property of said biochemical reaction network 
at said first time; 

(f) modifying said value provided to said variable constraint, and 

(g) repeating step (e), thereby determining a systemic property of said 
biochemical reaction network at a second time. 

72. The method of claim 71, wherein said value is modified based on 
said flux distribution at said first time. 

73. The method of claim 71, wherein said value is modified based on a 
change in an environmental condition. 
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74. The method of claim 71, further comprising repeating steps (e) 
through (g) for a specified number of timepoints. 
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